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We quantitatively study the transport of E. coli near the walls of confined microfluidic channels, 
and in more detail along the edges formed by the interception of two perpendicular walls. Our 
experiments establish the connection between bacteria motion at the flat surface and at the edges 
and demonstrate the robustness of the upstream motion at the edges. Upstream migration of E. 
coli at the edges is possible at much larger flow rates compared to motion at the flat surfaces. 
Interestingly, the bacteria speed at the edges mainly results from collisions between bacteria moving 
along this single line. We show that upstream motion not only takes place at the edge but also in an 
“edge boundary layer” whose size varies with the applied flow rate. We quantify the bacteria fluxes 
along the bottom walls and the edges and show that they result from both the transport velocity of 
bacteria and the decrease of surface concentration with increasing flow rate due to erosion processes. 
We rationalize our findings as a function of the local variations of the shear rate in the rectangular 
channels and hydrodynamic attractive forces between bacteria and walls. 


I. INTRODUCTION 

The question of motility and transfer of microorgan¬ 
isms in their environment is at the center of numerous 
issues. In many practical situations as in porous or frac¬ 
tured media, the pore size or the gap left for the micro¬ 
organisms to move can be very small and in theses con¬ 
fined situations surfaces become predominant. In addi¬ 
tion, internal walls of porous media, as well as biologi¬ 
cal conducts like blood vessels, lymphatic ducts, urinary 
and reproductive tracks, are not simply perfect cylin¬ 
ders, but their surfaces have irregularities, as grooves 
and crevices, that make them different from simple regu¬ 
lar surfaces. The understanding of bacteria motion along 
complex surfaces in the presence of flows is thus of strong 
importance for the control of microorganism transport in 
underground water resources, catheters or biological con¬ 
ducts. 

Bacteria transport has been investigated in bulk flows 
and Marcos et al. 0 have recently shown that bacteria 
drift with respect to the shear plane in Poiseuille flows. 
The drift has been explained as a consequence of the 
interaction between the chiral flagella of the bacteria and 
the shear rate 00. As a consequence bacteria drift 
in opposite directions in the upper and lower half of a 
Poiseuille flow, induced by the opposite signs of the shear 
rate components. 

The presence of walls modifies bacteria motion even 
without flow. In confined environments, bacteria are 
known to be attracted by flat surfaces and several stud¬ 
ies have shown that the concentration of bacteria is sig¬ 
nificantly larger at the top and bottom walls of square 
microchannels compared to the concentration in the bulk 
00 - 



FIG. 1. Sketch of the experimental configuration. A rect¬ 
angular microfluidic channel of height h = 20 pm and width 
w = 200 pm is observed using an inverted microscope. A 
bacteria suspension is flown at a volumetric rate Q in the di¬ 
rection indicated by the blue arrows. The red arrow in the x 
direction is our positive reference. Panels (a) and (b) indicate 
schematically typical bacteria motion at the top and bottom 
surfaces, respectively, as well as at their interceptions with 
the lateral edges. Edges with a green (red) color picture cor¬ 
respond to mostly “in-going” (“out-going”) mean transverse 
flux. 


The attraction towards the walls results from hydrody- 
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namic interactions between swimming bacteria and walls. 
Typically an autonomously swimming bacterium can be 
seen as a force dipole. Then, bacteria described as ’’push¬ 
ers” |1} are attracted to their specular hydrodynamic im¬ 
age close to a solid wall Furthermore, near a 

surface the rate of tumble has been observed to decrease 
with respect to the bulk @ also contributing to the long 
time bacteria spend very close to solid surfaces @,0. 
Bacteria motion at the surface is also modified compared 
to motion in the bulk due to lubrication forces between 
bacteria and walls. The viscous drag felt by the bacteria 
slows down the bacteria velocity 0 and leads to the ex¬ 
istence of circular trajectories due to their body rotation 
000. Note that a purely kinetic approach 00, 
not taking hydrodynamic interactions with walls into ac¬ 
count, also predicts increased bacteria concentrations at 
walls in confined geometries. 

Under flow, the interaction between local shear and 
swimmer motion close to the surface results in upstream 
migration at small shear rates and downstream swim¬ 
ming at a given angle with respect to the flow direction 
at higher shear rates, reported for E. coli in [15] and for 
mammalian sperm in [lfj. More generally, recent stud¬ 
ies suggest that upstream swimming takes place above a 
given threshold in shear rate for any front-back asymmet¬ 
ric micro-swimmer interacting hydrodynamically with a 
surface 0 - 

Concentration profiles stay flat in the bulk with a 
strong increase of concentration at the surfaces for small 
applied shear rates [18], but for higher shear rates more 
complex concentration profiles are observed and pre¬ 
dicted in the direction of the channel height [gj, 0, 0] . 
Ezhilan and Saintillan |0 and Chilukuri et al. 0 pre¬ 
dict a decrease of the surface concentration with increas¬ 
ing shear rate, but no experimental investigation of this 
phenomenon has been performed so far. 

Less work has been devoted to the study of bacte¬ 
ria motion at wall interceptions. Bacteria concentration 
have been found to be even higher at these edges com¬ 
pared to flat surfaces 0] and under flow bacteria motion 
is observed at the edges in a predominantly upstream 
way 0 , 0 ] over long distances. This upstream motion 
is at the origin of anomalous reconcentrations observed 
to be closely linked to the specific details of the confining 
structure 0,0. 

Here we study the behavior of bacteria while swim¬ 
ming in the vicinities and along the edges resulting from 
the interception between bottom and lateral walls of con¬ 
fined rectangular microchannels, as a first approach to 
understand their behavior in response to shear in con¬ 
fined irregular structures. 

We systematically quantify the concentrations and ve¬ 
locities of bacteria at the horizontal surfaces and at the 
edges, as a function of the applied mean shear rate. 
We show that concentrations decrease exponentially with 
shear rate due to erosion of bacteria. The slower decay 
of bacteria concentrations at the edges can be explained 
by the smaller local shear rate experienced by bacteria 


at the edges compared to the surfaces, as well as strong 
attractive hydrodynamic interactions and suppression of 
bacteria tumbles. We establish the link between bacte¬ 
ria motion at the horizontal surfaces and the anisotropic 
bacteria concentrations observed at opposite edges. We 
show that the bacteria navigation along the edges takes 
place at speeds mainly given by collisions between bac¬ 
teria moving along a single line along the edges. In con¬ 
trast, at the bottom and top surfaces bacteria are trans¬ 
ported downstream at a speed proportional to the mean 
shear rate, as soon as a critical shear rate is overcome. 

We define and measure an order parameter that allows 
to quantify the strong tendency of E. coli to navigate up¬ 
stream along the edges. It shows a transition from a sym¬ 
metric mix of downstream/upstream navigation at very 
low mean shear rates, to pure upstream navigation at 
larger mean shear rates. We identify an “edge boundary 
layer” close to the edges, where bacteria also navigate 
upstream. 

Finally, we show that the overall bacteria transport 
results from the dependency of both the concentrations 
and transport velocities on the mean shear rate. 


II. EXPERIMENTAL SET-UP AND METHODS 

Bacteria suspension - The bacteria are wild-type E. 
coli (ATCC 11105). Suspensions are prepared using 
the following protocol: 10 pL of bacteria were grown 
overnight in 15 mL of a rich culture medium (LB). 
From this, 5 mL was washed using 25 mL of PBS and 
centrifuged. Thereafter, the pellet was re-suspended 
into Minimal Motility Medium (MMA)0| and supple¬ 
mented with K-acetate (0.34 mM) and polyvinyl pyroli- 
done (PVP: 0.005%). After incubating for an hour in this 
medium to obtain a maximal activity, it was mixed with 
Percoll (1:1) to avoid bacteria sedimentation. This con¬ 
trolled environment promotes motility but does not allow 
bacteria replication. The experiments where performed 
at (25±2) °C. Under these conditions, the average swim¬ 
ming speed (far from the surfaces) is 28.7 ± 1.3 pm/s. 
E.coli perform run an tumble motions which can be re¬ 
garded as a directed random walk process 0, 0] . At 
short times, the trajectory is ballistic with a swimming 
velocity vq and at long times, the motion is diffusive with 
a translational diffusion coefficient: D = \v^t 1 where r 
is the cross-over time scale expressing the loss of direc¬ 
tional orientation and is r = 0.5s for our experiments 
0 , 0 . 

For the present experiments, the average bulk bacteria 
concentration is kept at = (3 ± 0.5) x 10 9 bact./mL 
(3 x 10~ 3 bact ./pm 3 ) unless otherwise stated. The 
volume fraction </>, based on a body volume of lpm 3 is 
cj) = 0.003 and corresponds to a dilute regime. 

MicroBuidic channel - The experimental cell is a 
rectangular channel made in PDMS using soft lithog¬ 
raphy techniques. The channel is h = 20 pm deep, 
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FIG. 2. Flow velocities and shear rates for a rectangular channel (w/h = 10) from the solution of Stoke’s equations (eqn. [TJ. 
(a) Normalized velocity field V(y, z)/Vmax- (b) Zoom of normalized velocity profiles close to the lateral edges at different z 
positions versus distance normalized by the channel half-width w/2. (c) Normalized shear rate amplitude 'yiocai{y,z)/ , y. (d) 
Zoom of the normalized shear rate profiles close to the lateral edges at different z positions versus distance normalized by the 
channel width w. 


w = 200 pm wide and several millimeters long (see 
Fig. a. Confinement of bacteria suspensions by two 
parallel walls can be quantified by comparing the 
distance between the two walls (the channel depth) to a 
typical distance over which bacteria swim between two 
successive tumbles, yielding the confinement parameter 
C = [H] . For C > la bacterium crosses the channel 

depth essentially without changing direction and the 
walls play an important role in the bacteria dynamics. 
On the other hand, for C « 1 the motion of bacteria 
within the channel depth is diffusive and the channel 
walls only play a minor role. For our experimental 
system, we find C = 0.75 and we are thus in a situation 
where walls play a non-negligible role. 

Experimental protocol and bacteria detection - The 
suspension is observed using an inverted microscope 
(Zeiss-Observer, Zl) with a high magnification objec¬ 
tive 100x (field-depth 6 pm). Videos are taken with 
a digital camera Photron FastCam SA3 (1024 x 1024 
pixels or 158 x 158 pm) at a frame-rate of 500 fps unless 
otherwise stated. The bacteria suspensions were seeded 
with latex beads (d = 2 pm, Beckman Coulter, density 
p = 1.027g/mL at a volume fraction 10 -5 %) and in¬ 
jected inside the micro-channel at different flow rates Q 
obtained by gravity over-pressure. The focus position is 
set at a height immediately above the bottom surface 
where a first homogeneous layer of moving bacteria can 
be detected. At this position the cell bodies appear as 


white areas surrounded by a dark halo. Through image 
post-processing they are detected via their local inten¬ 
sity maximum. Then using a calibrated criterion for the 
maximal intensity, we quantified that bacteria detected 
in this way where within a distance zq = 1.5 pm from the 
surface. Furthermore, we define that bacteria belong to 
the horizontal surfaces when they can be tracked in fo¬ 
cus for at least 1 second, which corresponds to a traveling 
distance without flow of typically their own size including 
the flagella bundle. 

Although the field depth of the microscope is quite 
small (6 pm), beads can be observed over the whole chan¬ 
nel depth when flowing through the micro-channel. The 
fastest bead in each video is used to determine the maxi¬ 
mum flow velocity V ma x in the channel. The uncertainty 
of this method depends on the number of observed beads 
and is almost negligible as soon as we observe more than 
10 beads. It can however be significant at very low flow 
rates. Maximal velocities detected vary between 0 and 
4150 pm/s. 


III. VELOCITY AND SHEAR PROFILES 

The rectangular channel we use here has an aspect ra¬ 
tio w/h = 10. Velocity profiles differ thus from an ideal 
parabolic Hele-Shaw flow profile, where the velocity is in¬ 
variant in the transverse direction y, for distances to the 
lateral wall larger than the cell height h. This will in- 
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fiuence the local values of flow velocities and shear rates 
and thus the transport properties of bacteria at the sur¬ 
faces. We present in Fig. [2] the flow velocities and shear 
rates computed for a rectangular channel with the same 
aspect ratio as our experimental channel. Velocities and 
shear-rates are normalised respectively by the maximal 
velocity V ma x and by the average shear rate 7 = 2V ™ a:L . 

The amplitude of the local shear rate is defined as: 
7 iocai(y,z) = + (^(y,z)) 2 . General views 

and zooms of the longitudinal velocities and shear rates 
are displayed on Figs. (2|a)(c) and Figs. [H(b)(d) respec¬ 
tively. The zooms show the regions close to the edges 
that we are specifically interested in. The flow and shear 
rate profiles are obtained from an exact solution of the 
Stokes equation for the longitudinal velocity field V(y, z) 
with no-slip boundary conditions [26| : 


V{y,z) = Y 

n,odd 
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where c = (with p the dynamic viscosity and p the 

pressure), — ^ < y < ^ and 0 < z < h. Interestingly, 
one can see that close to the edges there is a significant 
variation of the hydrodynamic conditions compared to 
the situation at the bottom/top surfaces. At a typical 
distance from the bottom wall, corresponding to our ob¬ 
servation plane ( z ~ 1.5 pm), we find a velocity reduction 
(Ysurface ~ 0.3 V ma x) and an increase in the local shear 
rate ( 7 surface ~ I. 77 ), while at the same height and lpm 
from the lateral wall, corresponding to the edge, we find 

Vedge ~ 0.0 5Vrnax and 7 edge ~ 0.57- 
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IV. RESULTS AND DISCUSSION 
A. Experimental observations 

In absence of flow, a large amount of bacteria swim 
close to the bottom surface and we observe them per¬ 
forming circular trajectories (no t shown) in agreement 
with previous observations [J, [ill, 113 • 

When a flow is applied, the scenario changes radically. 
The swimming direction is modified and at very small 
shear rates we observe upstream motion of bacteria, as 
also reported by Kaya and Koser [l5|- At higher shear 
rates (corresponding to most of our experiments) bac¬ 
teria navigate diagonally downstream (left or right as a 
function of the sign of the shear gradient), see Fig. [Hand 
snapshots on Fig. [3] This effect of “surface rheotaxis” is 
also in agreement with previous observations in jTB}. 

As bacteria reach an interception between the top and 
bottom surfaces with the vertical walls, they reorient 
along the edges, leading predominantly to an upstream 
motion. Due to shear, bacteria can also be eroded from 
the edges. In Fig. [H we present two typical trajectories 
corresponding to attachment (left panel) and detachment 


FIG. 3. Surface rheotaxis. Superposition of snapshots for the 
top and bottom channel surfaces. Individual bacteria trajec¬ 
tories in the channel are visible from image superposition and 
white arrows indicate the average direction of bacteria mo¬ 
tion. The green and red lines indicate respectively “in-going” 
and “out-going” edges, analogous to the ones sketched in Fig. 
[T] For these snapshots the bulk bacteria concentration of the 
reservoir was (4±0.5) x 10 9 bact./mL. The videos were taken 
at 40 x magnification and 30 fps with a digital camera Pix- 
eLINK PL-A741-E (512 x 512 pixels or 197 x 197 pm), at a 
mean shear rate of 19s -1 . 


(right panel) from an edge. Note that for these trajecto¬ 
ries bacteria swim indeed upstream at the edge. 

It is important to notice that surface rheotaxis has a 
strong influence on the balance of fluxes concerning the 
bacteria traffic at the 4 edges, as it breaks their geometri¬ 
cal symmetry. This is why we make a distinction between 
edges corresponding to “in-going” (green line) and “out¬ 
going” (red line), represented in Figs. [0 and [3] and that 
will be discussed further in section HVBI 
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FIG. 4. Attachment/detachment processes visualized at the 
edge. The background (25 x 50/rra) snapshots are taken near 
one lateral edge between the bottom surface and a lateral 
wall. The white line trajectories superimposed on the picture 
corresponds to bacteria detaching from the edge (left) and 
attaching to it (right). Mean shear rate is j = 220s _1 , (the 
edge corresponds to an “in-going” case, but it was not colored 
in green for clarity). 


B. Surface vs. edge erosion 


a 


Surface concentration 




Here we quantify the evolution of the bacteria concen¬ 
tration as a function of the mean shear rate at the flat 
horizontal surfaces and the edges. Fig. [5]shows the mean 
bacteria surface concentration at the bottom wall and at 
the edges of the channel as a function of the mean shear 
rate. 

In Fig. [5] the surface concentration a is defined as the 
number of bacteria observed within an area S of the sur¬ 
face, and within a distance Zo = 1.5/zto from it, divided 
by S. The volume concentration n s near the surface is 
then given by the number of bacteria observed in the 
same region, but divided by the volume zoS. The linear 
edge concentration A is defined as the number of bacte¬ 
ria observed along a distance L of the edge, divided by 
L. Finally, the volume concentration near the edge n e is 
given by the edge linear concentration divided by 1.5 [im 2 , 
i.e., we are assuming that bacteria swim along a corridor 
of width 1 fj.m and height 1.5 fim parallel to the edge. At 
zero shear rate, the volume concentrations for the surface 
and the edges shown in Fig. [5ja) and (b) are found to be 
around 1.5 xl0~ 2 and 10 _1 bact /\vrv? respectively. Com¬ 
pared to the “bulk” volume concentration of the reservoir 
of (rib = 3 x 10 -4 bact//zm 3 ) the concentration is thus 5 
times higher at the flat surfaces and 30 times higher at 
the edges. Ezhilan and Saintillan [l3|] predict a typical 
increase of the surface concentration at a flat surface for 
parameters comparable to our experimental conditions 
of typically 2-3 times. A purely kinetic model as theirs 
might thus not be enough to explain the large increase 
in concentration observed in our experiments. This indi¬ 
cates the importance of hydrodynamic attractive forces 


FIG. 5. Surface erosion: bacteria concentrations vs ab¬ 
solute shear rate value | 7 | . Lines are exponential fits: 
A exp (—I 7 I/ 7 0 ). (a) Bacteria concentrations are measured 

at the channel bottom surface, excluding 10 /rm-wide stripes 
from the lateral walls. A s = (2.81 ± 0 . 12 ) 10 _ 2 &act//rm 2 and 
7 ° = (143 ± 20)s _1 . (b) Linear concentrations at the edges 
for “in-going” and “out-going” edges . Fit parameters: A; = 
(2.0 ± 0.7)10 _ 1 bact//rm and j® = (250 ± 75)s _1 (in-going) 
and A 0 = (1.2 ± 0.7)10 ~ 1 bact/ and jo = (240 ± lOO )^^ 1 
(out-going). 


Q (or others) to maintain bacteria near the walls for long 
times. 

We now discuss the evolution of concentrations when 
a flow is applied. From Fig. [Sja) one can see that with 
increasing shear rate, the surface concentration decreases 
strongly, indicating erosion of bacteria from the wall. The 
concentration decrease is well described by an exponen¬ 
tial decay as A exp (—I7I/70) with 70 = (143 ± 20)s _1 . 
This typical erosion shear rate is large compared to the 
inverse of a typical hydrodynamic time scale for swim¬ 
ming bacteria, that can be defined as vq/1 , where vo is 
the bacteria swimming speed and l is a typical bacterium 
size. This difference might be explained once again by 
hydrodynamic interactions between the bacteria and the 
surface. 

The viscous drag on a bacterium can be estimated as 
F e oc with /i the fluid viscosity and where modifica¬ 
tions of the viscous drag due to the presence of the wall 
are neglected. On the other hand, the hydrodynamic at- 
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traction force in the presence of a flat surface scales as 
F a oc iil(—) 2 v o, with r w the distance between the swim¬ 
mer and the wall and l the dipole size, corresponding 
roughly to the bacteria body size (possibly including the 
flagella) Q. From the balance between F e and F a , a 
typical erosion shear rate can be obtained: 

7 ° oc v 0 /l. ( 2 ) 

Therefore, a distance of zo = 1 ym, l = 10/rm and 
v 0 = 30/zro/s yields: 7 0 « 3 x 10 2 s -1 , which is about 
the right order of magnitude for erosion, as seen in Fig. 
[5ja). However, this estimation is critically sensitive to 
the choice of the ratio — and clearly, should be vali- 
dated on a more refined hydrodynamic model, possibly 
also including the kinetic contribution of bacteria incom¬ 
ing from the bulk flow. 

Fig. 03(b) shows the linear bacteria density at the two 
edges corresponding to the bottom of the channel as a 
function of the mean shear rate. An important observa¬ 
tion can be made from this figure. The concentrations at 
the right and left edges are identical only at very small 
or very high shear rates. With increasing shear rate the 
concentration at the in-going edge increases whereas the 
concentration at the out-going edge decreases, leading to 
an asymmetry between the two edges. This has already 
been qualitatively observed in [271 and results from the 
bacteria transport at the flat surfaces. From the snap¬ 
shots on Fig. 03 it is clear that bacteria swim, in average, 
with a finite angle compared to the flow direction. This 
brings bacteria preferentially towards a given edge (the 
in-going edge). At the bottom surface bacteria are ob¬ 
served to drift towards the left with respect to the flow 
direction, leading to a concentration increase at the left 
edge, from the point of view of an observer moving with 
the flow. At higher shear rates a decrease of the bacte¬ 
ria concentration with shear rate is again observed. At 
very large shear rates the erosion (detachment) of bac¬ 
teria is so strong that the concentration tends towards 
zero for both edges. We attempt to adjust this decrease 
also by an exponential, and even if the quality of the fit 
for the edges is less good compared to the flat surface it 
leads to an estimate of the decay rate, which is found, by 
separately fitting the two data sets, to be approximately 
7 ° = (250 ± 100)s _1 for both edges. 

So, the shear rate associated with the concentration 
decrease is larger for the edges than for the surfaces and 
bacteria are thus eroded more slowly from the edges com¬ 
pared to the surfaces. In section Hill we have shown that 
the local shear rate at the edges is 3 — 4 times smaller 
compared to the local shear rate at the flat surfaces. This 
difference in local shear rate is enough to account for the 
difference in 70 observed. We can however, within our 
experimental resolution, not exclude other effects that 
might make bacteria more resistant to erosion at the 
edges, as increased hydrodynamic attraction at the cor¬ 
ner compared to the surface. Actually, the attractive in¬ 
teraction with the walls at the edges can be seen as stem¬ 


ming from the interaction between the bacterium and two 
specular images, each of them situated at the opposite 
side of each wall, plus a third specular image in a direc¬ 
tion extrapolated from the segment going from the actual 
bacteria to the interception between the walls. Further¬ 
more a smaller bacteria concentration at the flat surface 
for larger shear rates could also result in a smaller con¬ 
centration at the edges, as less bacteria reach the edges 
reducing bacteria concentrations there. 


C. Transport of bacteria by the flow 

1. At the flat surfaces 

We start by reporting the bacteria velocities at the 
flat surfaces. First, without flow, we noticed a signifi¬ 
cant decrease of bacteria velocities close to the surfaces 
(13 ± 2 fj,m/s) when compared with the velocities in the 
bulk (28.7±1.3 ym/s). Second, we report the transversal 
and longitudinal velocities of bacteria at the flat surfaces 
under flow. Bacteria move at the surfaces of the mi- 
crochannel following transversal, straight trajectories as 
is illustrated in Fig. [3] By tracking individual bacteria, 
we were able to obtain their velocities and orientation. 
We project the velocities on the x and y axes to obtain 
the longitudinal and transversal bacteria velocities. The 
transversal velocity on the surface always points in the 
vorticity direction, indicating that bacteria move to the 
left with respect to the flow direction at the bottom sur¬ 
face. The longitudinal direction at the surface is referred 
to the fluid direction, being positive when bacteria move 
downstream and negative when they migrate upstream. 

Fig. [bKa) shows the mean transversal bacteria velocity, 
oriented from the out-going to the in-going edge. If the 
flow is reversed, the transversal velocity also reverses, and 
the in-going and out-going edges exchange positions. The 
different symbols correspond to different flow directions. 
As the shear rate increases, we see a steep increase in 
the transversal velocity, until it saturates at a value that 
coincides with the average swimming velocity of the bac¬ 
teria at surfaces without flow. This means that bacteria 
swim almost perpendicular to the flow, as reported by 
Kaya and Koser [15|. This saturation occurs for a mean 
shear rate of the order of 40 s -1 . 

The bacteria moving at the surfaces feel a Stokes drag 
from the local flow and can thus be transported down¬ 
stream. In Fig. Ob), the mean longitudinal velocity far 
from the lateral walls is displayed as a function of the 
mean shear rate 7 . The effect of the flow is an entrain¬ 
ment proportional to 7 . Due to the fact that bacteria 
swim mostly in a direction perpendicular to the flow di¬ 
rection at higher mean shear rates the swimming speed 
does not influence the bacteria transport. The longitudi¬ 
nal bacteria velocity is found to be smaller than the local 
flow velocity {vi ongit udinal = 0 . 2 V sur / ace ), which can be 
explained by hydrodynamic interactions such as lubrica¬ 
tion forces slowing down the mean transport velocity. At 
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very small shear rates bacteria can swim upstream and 
average bacteria velocities in a direction opposite to the 
flow direction have been reported by Kaya and Koser 
pj|. Here we explore much higher flow rates, and the 
mean entrainment direction is essentially along the flow. 
Only for the smallest applied flow rate (see inset (S^b)) 
negative bacteria velocities has been detected. 
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FIG. 6. (a) Average transverse bacteria velocities at the hor¬ 
izontal surfaces as a function of the mean shear rate. Ex¬ 
perimental data collected for both positive (filled circles) and 
negative (open triangles) flow directions at the same surface. 
The horizontal line indicates the average swimming velocity 
of bacteria at the surface at zero shear rate, (b) Absolute 
value of the mean longitudinal bacteria velocities for various 
flow rates |'y | at the bottom horizontal surface. The line is 
a linear fit to the data with a slope A = 0.64 ± 0.02/rm. In 
inset, zoom on the longitudinal velocity < v x >. 


2. Upstream vs. downstream traffic at the edges 

In this section, we quantify in detail the traffic of E. 
coli moving along the edges. Fig[7Ka) displays the concen¬ 
tration of bacteria swimming upstream and downstream 
at a given edge as a function of the mean shear rate. In¬ 
terestingly, the number of bacteria swimming upstream 
along the edges largely dominates the number of bac¬ 
teria swimming downstream. The slight asymmetry in 
the curves for positive and negative mean shear rates is 
due to the fact that the edge switches from an ingoing 
edge (7 < 0 ) to an outgoing (7 > 0 ) edge when the flow is 
inverted. The concentration of bacteria swimming down¬ 
stream is observed to be very small as soon as the mean 
shear rate is larger than |^y | > 25 s -1 . Shear induced 


orientation near the edges seems to be the cause of this 
difference: as bacteria approach the edge transversally as 
illustrated in Figs. [I] and 01 the local shear rotates the 
bacteria body that aligns preferentially in the direction 
facing the flow [ 2 ^, [2lJ . This shear induced orientation of 
bacteria does not take place at very small flow rates, and 
the concentrations of upstream and downstream bacteria 
at the edges are similar. 

In order to further quantify the bias in the direction of 
the bacteria traffic along the edges, we define a parameter 
called Bias in the Edge Traffic or BET number, as: 


BET = 


A up A down 
A up T A down 


(3) 


where X up and A down are respectively the linear concen¬ 
trations of bacteria moving in the positive and negative 
directions relative to the red arrow shown in Figs. [I] and 
01(a). BET thus allows to visualize the bias in the bac¬ 
teria navigation without taking the total bacteria con¬ 
centration or the erosion of bacteria from the edges into 
account. Fig. (lib) shows the dependence of BET on the 
mean shear rate 7 for different experiments, i.e. differ¬ 
ent channels of similar geometry, different concentrations 
and measured at different edges. When 7 = 0, there is 
an equal amount of bacteria moving up and downstream, 
so BET = 0. As the shear rate increases in the positive 
direction, more and more bacteria move upstream (i.e., 
in the negative direction), until BET = — 1. A similar 
reasoning explains the shape of the curve of Fig. [7Kb) for 
negative shear rates. 

We now consider bacteria velocities at the edges. The 
average velocities of bacteria moving upstream and down¬ 
stream at the two bottom edges are displayed on Fig. 
□(c). Bacteria velocities are found to be identical for the 
two populations for very small shear rates. Velocities of 
about lOfim/s (slightly smaller than the swimming speed 
at the flat surface) are observed. Surprisingly, with in¬ 
creasing shear rate upstream swimming bacteria become 
faster in average, whereas downstream swimming bacte¬ 
ria slow down. At even higher shear rates, the tendency 
is reversed and the downstream swimming bacteria see 
their velocity increased whereas the upstream swimming 
bacteria see their velocity decreased. These observations 
indicate transport dynamics very different from those re¬ 
ported at the flat surface, where the mean bacteria ve¬ 
locity is directly proportional to the shear rate (see. Fig. 
[ 6 Kb)). The transport velocities at the edges can thus not 
be simply explained by the hydrodynamic drag exerted 
by the local flow, but we associate these observations 
with the specific dynamics of bacteria moving along a 
unidimensional corridor along the edges. 

In general, bacteria can move up and downstream at 
the edge, leading to collisions between swimmers. We 
have observed that during frontal collisions, bacteria slow 
down for a certain time until they cross each other and 
continue swimming along the edge. We have also ob¬ 
served bacteria to detach from the edge during such a 
collision. Similar slowing down or detachment can also 












a 



Y (s' 1 ) 


b 



c 



IyI (s' 1 ) 


FIG. 7. Bacteria transport at the edges, (a) Concentration 
of bacteria along the left edge. Note that the “in-going” edge 
for 7 mean < 0 becomes an “out-going” edge when 7 meQ „ > 0 . 
(c) Bias parameter in Edge Traffic (BET) as a function of 
the shear rates (see Eq. [3]), The meaning of the labels in the 
different symbols in the picture are: LI and R1 Left and right 
walls respectively in an experiment at a bulk concentration 
in the reservoir of (3 ± 0.5) x 10 s bact./mL; L2 left wall at 
(4 ± 0.5) x 10 s bact./mL. (c) Absolute value of the mean 
longitudinal velocity for the bacteria population moving along 
the edges. The horizontal line is their mean velocity at the 
center of the channel at zero flow. 


take place during rear collisions between two bacteria 
swimming in the same direction, but are much less fre¬ 
quent. 

For very small shear rates the concentration of bacteria 
at the edges is large, and collisions are frequent. These 
collisions are at the origin of the decreased transport 


speed of bacteria at the edges compared to the flat sur¬ 
faces. As soon as the shear rate is larger than 25s -1 the 
number of bacteria swimming upstream is much larger 
than the number of downstream swimmers (Fig. Eta)), 
resulting in more frequent collisions for bacteria swim¬ 
ming in the flow direction, and the downstream swim¬ 
ming population is thus more strongly slowed down, pro¬ 
voking a decrease of their velocity compared to the up¬ 
stream swimming population. 

With increasing flow rate, the concentration of bacte¬ 
ria at the edges decreases and bacteria transport veloc¬ 
ities increase again, due to a decrease in the number of 
collisions. At very high shear rates, only very little bac¬ 
teria are left at the edges and we expect collisions not to 
play an important role any more. In this range of mean 
shear rates, the mean upstream bacteria velocities are 
observed to be close to the swimming speed measured 
at the flat surfaces without flow, indicating that bacteria 
swim upstream nearly undisturbed by the local flow. It 
is worth noticing that the speeds of individual upstream 
and downstream bacteria between collisions are always 
roughly similar and also very close to the bacteria swim¬ 
ming speed measured at the flat surfaces without flow. 
For the downstream swimming bacteria the small aver¬ 
age velocities measured in this range of shear rates are 
due to their attachment dynamics. During a typical at¬ 
tachment process under flow, a bacterium, advected by 
the main stream arrives at an edge, flips and starts swim¬ 
ming upstream. During this flipping process, the swim¬ 
mer changes its velocity direction from downstream to 
upstream, contributing to the statistics of downstream 
swimmers velocity for a short lapse of time, with a veloc¬ 
ity that decreases to zero. 

At even higher shear rates the viscous drag on the bac¬ 
teria might become important, leading to a decrease of 
the mean velocity of upstream swimming bacteria and 
a decrease for the downstream swimming bacteria. At 
a shear rate of 7=400s _1 the local flow velocity at the 
edge is 20/im/s. At the flat surfaces bacteria where trans¬ 
ported downstream at velocities five times smaller com¬ 
pared to the local flow velocity. Assuming a similar re¬ 
lation at the edges this would correspond to a transport 
velocity of 4 fim/s that needs to be added to the swim¬ 
ming speed of the bacteria. This is not in contradiction 
with the observed decrease in average transport veloc¬ 
ity of upstream swimming bacteria at these high mean 
shear rates. Most likely a correct modeling of the lu¬ 
brication forces would lead to an even smaller expected 
transport velocity at the edges. The striking fact that 
bacteria move at the edges nearly unperturbed by the 
flow can thus at last partially be explained by the de¬ 
creased local flow velocities at the edges compared to the 
flat surfaces. Increased drag close to a corner and the 
significant strength of the hydrodynamic interactions be¬ 
tween bacteria and the edges might even enhance this 
effect. 

As overall bacteria speeds at the edges vary little as 
compared to concentrations variations, the total bacteria 
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flux at the edges is directly proportional to the concen¬ 
tration profile represented on Fig. [7[a) and has not been 
represented separately. 


D. Edge boundary layer 


The previous measurements show that the edges have 
singular transport properties, as there is a significant flux 
of bacteria moving against the flow. Here we characterize 
whether this upstream motion is restricted to the edges, 
by measuring the bacteria flux along the flow direction 
at the surfaces and close to the edges and identify what 
we call the Edge Boundary Layer (EBL). 

The bacteria flux at the bottom surface is represented 
on Fig. [Ha) as a function of the distance from the wall. 
If bacteria are transported upstream, the bacteria flux is 
negative, and it is positive if they are transported down¬ 
stream. As can be seen in Fig. (Ha), at relatively large 
distances from the edge, bacteria are advected down¬ 
stream with the flow and only at very small shear rates 
(see results for 7 = Is -1 ) they move upstream even far 
away from the edge. The bacteria flux far away from 
the edges, the saturation flux J 0 , is represented in Fig. 
(Hb). When increasing the flow rate, the bacteria flux 
first increases and then decreases again, illustrating the 
interplay between bacteria transport and erosion. First 
the flux increases due to the linear increase in flow ve¬ 
locity (Fig. (Hb)), then, the effect of erosion starts to 
dominate and exponentially decreases the surface concen¬ 
tration (Fig. [Ha)), leading to an overall decrease of the 
surface flux. This can be quantified comparing the mea¬ 
sured saturation flux Jo with the product a (vi ong it u dinai) 
from the expressions obtained by fitting the experimental 
data in Fig. (Ha) and Fig. (Hb). The result is shown in 
Fig. (Hb) as the black line. We can also adjust the re¬ 
sults for Jo vs the mean shear rate directly with the same 
functional dependence Jo = Ayexp — 7 / 70 - The value 
obtained for y c is within the error bars in agreement with 
the critical shear rate for erosion at the flat surfaces 7 ° 
from Fig. Ha). From the fit of the bacteria density and 
velocity we obtain for A = 1.8 x 10 - 1 bact//im slightly 
smaller than the prefactor obtained from the best fit to 
the saturation flux. 

Close to the edge bacteria are observed to move up¬ 
stream even at mean shear rates where the flux far from 
the edges is observed to be downstream. From the longi¬ 
tudinal flux profile, we determine the distance from the 
edge where the flux changes sign, defining the width of 
the EBL, which is shown in Fig. He). The boundary 
layer builds up when the flow is turned on and reaches 
a maximal width of about 10 /rm at a mean shear rate 
of about 7 = 20s -1 . When further increasing the mean 
shear rate the EBL decreases again and stabilizes at a 
value of 2fj,m corresponding to the width of a single bac¬ 
terium and bacteria move upstream in a single line along 
the edge. 


a 



b 



I Y I (s' 1 ) 


c 



I Y I (s' 1 ) 


FIG. 8 . Edge Boundary layer: (a) longitudinal flux density 
of bacteria J x as a function of the vertical direction y for 3 
values of the shear rate 7 . For negative values of the hori¬ 
zontal axis bacteria move upstream inside the edge boundary 
layer, while they move downstream for positive values of the 
horizontal axis. The continuous lines are exponential fits of 
the form Jo + Bexp (— y/5 ) where Jo, B and 5 depend on 7 . 

(b) Longitudinal Bacteria flux Jo from the exponential fitting 
of each bacteria flux vs 7 . The black solid line is the prod¬ 
uct a (viongitudinai ) from independent fittings and the blue 
dashed line is the best fit to the data using Ay exp — y/ 7 c , 
where A = (2.5 ± 0.3)10 - 2 bact//rm and y c = 110 ± 10s -1 . 

(c) Width of the Edge Boundary Layer as a function of the 
mean shear rate 7 . 
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V. CONCLUSIONS 

In this paper, we have quantified the transport of bac¬ 
teria at flat surfaces and edges in confined microfluidic 
channels under flow. 

We have measured the decrease of the surface and edge 
concentrations as a function of applied mean shear rate. 
The slower decrease of the bacteria concentration at the 
edge can be explained by the smaller local shear rate at 
the edge and possibly increased hydrodynamic interac¬ 
tions at the intersection between two flat surfaces. Fur¬ 
thermore bacteria concentrations at different edges are 
not identical, because rheotaxis at the horizontal surfaces 
breaks the symmetry along the main axis of the channel, 
bringing bacteria preferentially towards a given edge: for 
a given shear flow in one direction, two diagonally op¬ 
posed edges of the rectangular cross-section are outgoing 
edges, and the other two are ingoing edges. 

We have observed that bacteria swim predominantly 
upstream at the edges as soon as a small flow is applied. 
This is attributed to shear induced reorientation of the 
bacteria attaching to the lateral edges. We have quanti¬ 
fied the strong bias in the swimming direction of bacte¬ 
ria at the edges towards upstream swimming by means 
of an order parameter that accounts for the symmetry of 
up and downstream swimming bacteria concentrations at 
the edges: the proportion of bacteria moving upstream 
increases very quickly as the shear rate increases until all 
bacteria are observed to swim upstream. 

Bacteria swimming at the edges are nearly undisturbed 
by the applied flow even at mean shear rates where bac¬ 


teria transport at the flat surfaces is already strongly 
influenced by the latter. This can be explained again 
by the decreased local shear rate at the edges and pos¬ 
sibly increased hydrodynamic interactions. Interestingly, 
we have found that average bacteria velocities along the 
unidimensional corridor of the edges are mainly the re¬ 
sult of collisions between up and downstream swimming 
bacteria. The number of collisions depends on the total 
concentration of bacteria at the edges as well as the per¬ 
centage of up and downstream swimming bacteria and is 
thus found to be a non monotonic function of the shear 
rate. 

Bacteria are able to swim upstream not only at the 
edges, but as well within an edge boundary layer (EBL). 
The width of the EBL decreases from approximately 10 
bacteria body widths at small shear rates, to 1 body 
width at higher shear rates. 

Our results thus quantify the bacteria fluxes at all sur¬ 
faces of a confined microchannel as a function of the mean 
shear rate. In the future these results can be used to un¬ 
derstand bacteria transport in more complex geometries 
or to design specific flow geometries to guide bacteria 
fluxes to selected positions. 
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